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We study the order parameter probability distribution at the critical point 
for the three-dimensional spin-1/2 and spin-1 Ising models on the simple cu- 
bic lattice with periodic boundary conditions. The finite size scaling relation 
for the order parameter probability distribution is tested and verified numer- 
ically by microcanonical Creutz cellular automata simulations. The state 
critical exponent 5, which characteries the far tail regime of the scaling order 
parameter probability distribution, is estimated for 3-d Ising models using 
the cellular automaton simulations at the critical temperature. The results 
are in good agreement with the monte carlo calculations. 
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A quantity of central importance for the finite-size scaling analyses of the 
critical phenomena is the order parameter probability distribution P(M). 
The knowledge of the finite size scaling function for P(M) of the Ising 
model makes it possible to calculate all the moments of the order parameter 
and all its cumulants.' 1-5 ' Most properties of P(M) are known from computer 
simulations J 6 _14 1 For d=3, the corresponding finite size scaling function do 
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not exist in the analytical form. There are effort to get simple analytical 
functions by fitting to Monte Carlo results at the critical temperature for the 
infinite lattice. I 10-13 ! The main purpose of these study is to test the finite 
size scaling relations for the order parameter probability distribution of the 
three-dimensional spin-1/2 and spin-1 Ising models by the microcanonical 
Creutz cellular automaton algorithm on the simple cubic lattice and to obtain 
the value of the state critical exponent 8 . The Creutz cellular automaton 
(CCA) for the Ising model has been proven to be successful in producing the 
values of the universal statical critical exponents and the critical temperature 
in two and higher dimensions J 15-20 ] The CCA algorithm, which was first 
introduced by Creutz, I 20 ! is a microcanonical algorithm interpolating between 
the canonical Monte Carlo and Molecular dynamics techniques. 

In this paper, the probability distribution of order parameter is obtained 
for the two variants of the Ising model on the CCA. The first model is the 
3-d spin-1/2 Ising model on the simple cubic lattice. For the zero external 
magnetic field (H — 0), the Hamiltonian of the model is given by 

Hi = - J E S * S > 

<ij> 

where Si = ±1 and the sum is carried out over all nearest-neighboring (nn) 
spin pairs. The parameter J ( J > ) is the ferromagnetic coupling constant. 
The simulations are carried out on simple cubic lattice LxLxL of linear di- 
mensions L = 16, 18, 20, 24 and 40 with periodic boundary conditions. The 
second model is the Blume-Capel(BC) model without single-ion anisotropy 
parameter ( D — ). The BC model is a spin-1 Ising model.' 21 ' 22 ] It has 
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the same Hamiltonian with spin-1/2 Ising model for the H = . Here, the 
spins can take three discrete values -1, and 1. The BC model is simulated 
using an improved algorithm' 18 ' from CCA for simple cubic lattice LxLxL of 
linear dimensions L = 16, 20, 24, 28 and 32 with periodic boundary condi- 
tions. The data are averages over the lattice and the number of time steps 
(1.000.000) during which the cellular automaton develops. The simulations 
are done 20 times with different initial configurations at the critical point for 
the Ising models. 

The infinite lattice critical temperature values for spin-1/2 and spin-1 
Ising models are estimated from the temperature variation of the Binder 
forth-order cumulant for the finite lattices. The Binder forth-order cumulant 
of the order parameter, which is used to estimate the infinite lattice critical 
temperature, is given by 

9l-1- (2) 

The temperature variations of the Binder cumulant are illustrated in Fig.l 
(a) for spin-1/2 and (b) in Fig.l(b) for spin-1 Ising model. The infinite lattice 
critical temperatures are obtained from the intersection of the finite lattice 
Binder cumulant curves. The critical temperature value for spin-1/2 Ising 
model is estimated as kT c / J = 4.511±0.002 which is in good agreement with 
Monte Carlo calculations J 23-26 ' For the spin-1 Ising model , the estimated 
critical temperature (kT c / J = 3.197 ± 0.002) is also in a good agreement 
with a series expansion results. I 27 " 29 ! The behavior of the order parameter 
probability distribution at the critical point is investigated at these estimated 
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critical temperatures. 

The order parameter probability distribution Pl(M) is calculated by 



Pl{M) 



N m 



(3) 



Nccas 



where Nm is the number of times that magnetization M appears, and Nccas 
is total number of Creutz cellular automaton steps. Histograms of 200 bins 
are used in plotting P L (M) at critical point for finite lattice sizes which are 
shown in Fig.2 (a) for spin-1/2 and in Fig.2(b) for spin-1 Ising model. The 
simulations have been performed for the twenty different initial configura- 
tions with a constant total energy at critical point. The data of Pl(M) are 
obtained from averages of data for different initial configurations. The rela- 
tive standard deviation (RSD) of the average Pl(M) values is approximately 
in the interval 5% — 10%. 

The finite size scaling hypothesis for the probability distribution of order 
parameter for Ising model can be expressed generally as follows^ 



and P*(M*) is a universal scaling function. To compute the normalized 
distribution P*(M*) via Eq.4 one has evaluate the pre-factor b. The value 
of b for each lattice sizes can be easily calculated by 



P L (M) = bP*(M*) 



(4) 



where b = b L^ u , (3 and v are critical exponents, M* = bM, b is a constant, 




(5) 
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Thus, we used the Eq.4 for obtained of the universal function P*(M*) as 
in Ref. [8]. At the critical point, the plots of finite size scaling for P(M) 
are illustrated in Fig.3 (a) for spin-1/2 and in Fig. 3(b) for spin-1 models. 
The micro canonical simulations have been done on simple cubic lattices at 
kT c /J = 4.511 for spin-1/2 and at kT c /J = 3.197 for spin-1. For both mod- 
els, the universality at critical point can be easily seen from Fig. 3 (a) and (b). 
The scaling probability distributions P*(M*) for spin-1/2 and spin-1 Ising 
models coincided with each other for all M* values. These results verify the 
finite size scaling relation given in Eq. 4 for the order parameter probability 
distribution of the three dimensional Ising model. On the other hand, the 
log-log plots of the pre-factor b against L are shown in Fig. 4. The slopes 
of the data line a single curve are given the values of /3/z/=0.498±0.003 for 
spin-1/2 and /3/z/=0.492±0.003 for spin-1. Furthermore, the (3/v values for 
both models are estimated using the scaling relation of the order parame- 
ter at the critical temperature [Fig.4(b)]. The slopes of the log-log plots of 
M(T C ) against to L are given the values the /9/z/=0.50±0.01 for spin-1/2 and 
/3/z/=0.51±0.01 for spin-1. All estimated values of (3/v are in good agreement 
with universal valued' 30 !. 

On the other hand , for the \M*\ > 1 and the periodic boundary condi- 
tions, the scaling function of the order parameter probability distribution is 
expected to have the following exponential form^ 13 ' 

P*(M*) oc exp {-AM* 5+1 ) (6) 

where 5 denotes the equation of state exponent. Its value can be obtained 
from the log-log plot of lnP*(M*) against to M*. In the far tail regions, the 
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slopes of data line a single curve gives approximately 5 + 1 = 5.8 for spin-1/2 
and 5 + 1 = 5.7 for spin-1 on the right and left tails. The behavior of tails 
for large M* are shown in Fig. 5 (a) and (b) for spin models. The estimated 
values of 5 + 1 are in good agreement with Monte Carlo results (5 = 4.8) ^' 13 ^ 
for 3-d Ising model. Finally, the finite size scaling relation for the order 
parameter probability distributions of the Ising models is verified by Creutz 
cellular automaton simulations numerically. 
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Figure Captions 

Fig.l The variations of Binder cumulant against to kT / J for (a) spin-1/2 
and (b).spin-l Ising model. 

Fig. 2. The order parameter probability distribution at the critical point 
for (a) spin-1/2 and (b).spin-l Ising model on the simple cubic lattices. Simu- 
lations were performed at kT c / J = 4.511 for spin-1/2 and at kT c / J = 3.197 
for spin-1 model. 

Fig. 3. Scaling functions P*(M*) for (a) three-dimensional spin-1/2 and 
(b) spin-1 Ising model on simple cubic lattices. 

Fig. 4. (a ) The Log-log plot of b against to L for spin-1/2 and spin-1 Ising 
models, (b) The log-log plots of M(T C ) against to L . The values of slopes 
are in good agreement with the universal (3/v value. 

Fig.5. The log-log plots of ln( P*(M*) ) against to M* for (a) spin-1/2 
and (b) spin-1 Ising model. The slopes are equal to 5 + 1 = 5.8 for spin-1/2 
and 5 + 1 = 5.7 for spin-1 model in the far tail regions. 
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